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Abstract 

Wireless networks found in practice are finite, and hence non-stationary, with nodes typically 
non-homogeneously deployed over the network area. This obviously leads to a location-dependent 
performance and to boundary effects which is often neglected in network modeling. In this work, 
interference and local throughput in a Poisson network, where the spatial distribution of nodes is 
isotropic but not necessarily stationary, are studied. They are precisely analyzed as a function of (i) 
an arbitrary receiver location and of (ii) an arbitrary isotropic shape of the spatial distribution. The 
interference distribution is characterized through a first moment analysis for arbitrary block-fading 
channels (including the pure path loss model) and bounds on the tail probability are derived. For 
Rayleigh fading, the Laplace transform of the interference distribution is presented. For the path losses 
a = 2 and a = 4 closed-form results are derived. Two metrics suitable for measuring local throughput 
in non-stationary networks are proposed, and they are discussed for the isotropic model at hand. One 
the one hand, this work reveals some interesting and fundamental facts, particularly it revises some 
prior results for the case a = 2. On the other, it provides a powerful tool for studying non-stationary 
networks as demonstrated through numerous examples. 
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I. Introduction 

Stochastic Geometry, in particular the theory of point processes, has attracted much attention 
in the field of interference modeling and performance evaluation for wireless networks with 
(multiple) sources of randomness, such as user locations, user behavior and channel fading. 
Originally inspired by statistical problems in material sciences [1], it has enabled researchers 
to answer some of the fundamental questions in wireless network theory; thereby bridging the 
gap between yet unsolved multi-user communication problems in information theory and design 
guidelines for network operators. Among the numerous contributions in this area, the perhaps 
most notable one has been the transmission capacity framework, which led to many profound 
results in the topic of wireless networks, cf. ED-flTl. 

In essence, the transmission capacity framework models the spatial distribution of nodes of the 
network as being random rather than assuming a deterministic or fixed spatial configuration. The 
advantages hereof are manifold, see e.g. [[U, flU for further details; but perhaps most significantly, 
such a probabilistic approach decouples the network model from the actual spatial configuration 
which, in turn, increases the generality of obtainable results. Moreover, these results can be given 
in closed-form, thereby revealing the individual effects of the system parameters on the network 
performance. By definition, the transmission capacity gives the maximum density of successful 
transmissions that can take place simultaneously in the network, weighted by the probability of 
success of a typical transmission H|. 

With only a few exceptions the node locations have mostly been modeled by a stationary 
point process, i.e., the statistical properties^ do not depend on the particular location in the 
network. Stationarity is a desirable property, since it allows for analytic tractability and, even 
more important, it represents a key requirement for applying the transmission capacity metric. 

In practice though, wireless networks exhibit a non- stationary spatial node distribution; the 
perhaps most obvious example is when the network area is finite. In such a network the 
interference situation near the center significantly differs from that at the boundary. Besides 
this simple example, more complex deployments are often found in practice, e.g., wireless 
sensor networks created by airdrop [8] or spontaneous formation of hot spots [9|. The spatial 
configuration of such hotspots is typically dictated by user motion and by geometric constraints 

'Such as interference, outage probability or throughput. 
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as illustrated in [TTOll for the example of a campus-wide Wi-Fi deployment. There is furthermore 
a growing need for cellular operators to better understand not only the temporal variations in 
user traffic demands but also its spatial dependence; the optimal interplay between small-cell 
deployment and Wi-Fi off-loading -a promising approach to boost capacity in dense areas- 
requires carefully pinpointing areas of peak-traffic demands ifTTTl and identifying locations for 
deployment so to reach the mobile users |fT2|. In this context, analytic tools for quantifying the 
network performance while taking into account user mobility as well as hard-to-predict spatial 
configurations are of eminent importance. 

Even though the stationarity assumption has not really narrowed the range of obtainable 
insights, it has some shortcomings to the analysis of wireless networks which are described in 
the following: 

A. Limitations of the stationarity assumption 

1 ) Infinite networks: Stationarity implies that the network consists of infinitely many nodes 
spread over an infinitely large region. However, the number of nodes as well as the network 
area is finite in practice. This discrepancy for example prevents the use of absolute performance 
metrics, e.g., the sum capacity, since the number of transmissions is unbounded in this case. 

2) Boundary effects: Boundary effects are inherently ignored when assuming stationarity. 
Certainly, boundary effects play a critical role in real-world networks as they cause heterogeneity 
in the nodes' capabilities, i.e., nodes being dis-/connected, interference-moise-limited, etc. A 
careful treatment of boundary effects is crucial to network design in order to assure the level of 
quality of service (QoS) that was targeted before deployment. 

3) Model artifacts for path loss exponent a = 2: In the stationary case and with a path loss 
exponent of a = 2, the sum interference is almost surely (a.s.) infinite [3J. This model artifact 
results in an outage probability of one or, equivalently, to a transmission capacity equal to zero. 
More specifically, stationary models lose their accuracy as the path loss exponent tends to a = 2 
because the sum interference becomes dominated by the infinite number of far nodes. 

4) Application of the transmission capacity metric: As already mentioned, the transmission 
capacity metric only applies to stationary networks. The reason is that this metric makes a 
global statement about the network based on a local analysis. Hence when the local analysis 
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is location-dependent, it is no longer possible to infer the global performance from the local 
analysis. 

B. Contribution 

With this work, we extend the already existing framework by relaxing the stationarity assump- 
tion for the Poisson point process, thereby allowing for isotropy only. Note that we retain the 
Poisson property of the point process. The stationary case is included in our results as a special 
case. Specifically, we have the following contributions: 

• The interference statistics (first moment and distribution function) are analyzed for both 
arbitrary and Rayleigh block-fading channels. In the case of arbitrary fading, we derive in 
closed-form the exact first moment as a function of the spatial shape of the node distribution 
and of the arbitrary receiver location for the path loss exponents a = 2 and a — 4. 
Using these results, an upper bound on the tail probability of the interference is derived. A 
corresponding lower bound that is not limited to the above values of a is also presented 
for suitable spatial shapes. In case of Rayleigh fading, we derive the Laplace transform of 
the interference which can be used for computing outage probabilities, also in closed-form, 
as a function of the spatial shape and of the arbitrary receiver location for the above values 
of a. 

• We demonstrate how known results from the literature for stationary Poisson point processes 
arise from our results as a special case. One important insight is that for a — 2, one can 
find cases where a.s. infinitely many nodes contribute to the sum interference while the sum 
interference remains a.s. finite. This result sharpens prior statements about the nature of the 
interference for a — 2, and suggests that there exists a transition between sparse and dense 
networks. 

• We evaluate and demonstrate the use of our exact model by comparing it to a stationary 
model that uses a local approximation to capture the non-homogeneity in the spatial distribu- 
tion of nodes. We show that, depending on the spatial shape in question, large approximation 
errors are avoided by using our exact model, particularly in transition areas where boundary 
effects should be properly captured. 

• We further discuss how to measure throughput in non-stationary networks and propose 
two metrics that are capable of quantifying the local throughput in such networks while 
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taking into account their non-stationary nature. The first one is called local transmission 
capacity, which is a refinement of the original transmission capacity metric. It provides 
spatial information about the local throughput around a certain point. The second one is 
called sum spatial throughput and resembles the spatial throughput metric, which yields 
the density of successful transmissions in a stationary network. In contrast, the sum spatial 
throughput gives the absolute (sum) number of successful transmissions, thereby accounting 
for all non-homogeneities and boundary effects of the network. The two metrics proposed 
in this work are introduced for more general non- stationary networks and are furthermore 
discussed for the isotropic Poisson model at hand. 
• We give several examples and highlight some applications to which the developed model 
can be applied. In particular, we discuss how the accuracy of models for carrier-sense 
based decentralized networks can be significantly increased by means of our results. Also 
interestingly, for a — v 2 our analysis fundamentally extends prior work dealing with the 
famous question whether or not to use spread spectrum techniques for increasing network 
throughput. 



C. Related work 

Interference modeling and network analysis using tools from Stochastic Geometry has become 
a multi-faceted research field [|2]|-[|7]|, [fT3l - ll23l . These works mostly differ by the types of 
applications under consideration. They may furthermore be categorized by whether the underlying 
spatial model is stationary (homogeneous or non-homogeneous) or non- stationary. 

Early works dealing with interference modeling assumed a stationary^ Poisson point process 
for the interferer locations, cf. 11131 . |fl4l. The transmission capacity framework -and the many 
fundamental results that arose thereof- dealt initially with the same kind of point processes 
j3l-[|7l. These works enormously contributed to a better understanding of the basic interactions 
between system parameters of a wireless network. 

Stationary models with non-homogeneous node deployments have already been investigated. 
Specifically, Poisson-Cluster |[T5l and Matern hard-core models lfT6l were studied, as they are 
well-suited for analyzing more sophisticated medium access control (MAC) schemes. Treated as 

2 For the Poisson point process stationarity implies homogeneity, i.e., constant density [24|. 
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general motion-invariant, these and similar models were further analyzed in 0, IfTTl , |[T8l in a 
unifying way. Finite networks may be modeled by stationary cluster point processes and letting 
the parent density tend to zero (with the representative cluster left). However, the problem of 
infinite interference for a = 2 is not solved by that as with all stationary models. 

The need for non-stationary models for characterizing more complex node deployments in 
wireless networks was reported in |fT9l . The authors discussed techniques that generate non- 
uniform node distributions for the purpose of efficient network simulations. In [20J, a non- 
stationary and isotropic node distribution was assumed for analyzing multi-antenna receivers. 
While the analysis showed that the shape of the spatial distribution has a considerable impact on 
link performance, the scenario was limited only to the case of the receiver located in the origin. 
In ll25l . a first attempt was made towards analyzing the interference, resulting from an non- 
stationary Poisson point process, at an arbitrary receiver location and for an arbitrary isotropic 
node distribution. 

D. Notation and paper organization 

Sans- serif- style letters (z) represent random variables while serif-style letters (z) represent 
deterministic variables or constants. We denote by P(-) and E [•] the probability measure and 
the expectation operator, respectively. We will sometimes include the random variable under 
consideration in the subscript whenever necessary, e.g., P z (■) and E z [■]. The imaginary unit is 
given by j = \J—1 and £H(-) denotes the real part of a complex-valued number. The Euclidean 
distance is denoted by | • |. We define by b(z, r) and db(z,r) a disc and a circle, respectively, 
centered at z E M? with radius r > 0. The origin (0, 0) is denoted by o. 

The paper is organized as follows: Section [II] introduces the network model and provides a 
discussion about the assumptions made in this work. In Section |inj the interference statistics are 
analyzed as a function of an arbitrary location in the network and of an arbitrary isotropic spatial 
distribution of nodes. Here, we first focus on arbitrary fading channels, including the pure path 
loss case, and derive first moment results and bounds on the tail probability of the interference. 
Then, we treat the special case of Rayleigh fading and compute the Laplace transform of the 
interference distribution. These results are then used in Section[IV]to study the outage probability. 
Here, we also compare our exact model to a local approximation in order to quantify the gain in 
accuracy that we obtain by our model. Section |V] is devoted to the question of how to measure 
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local throughput in non- stationary networks. We give some discussions on that and show how 
the application of these metrics simplifies in the special case of isotropy. We also present some 
examples and applications to which our results can be gratefully applied. Section |VT] concludes 
the paper. 

II. Network Model and Assumptions 

We consider a packet-based wireless network with identically-equipped nodes isotropically 
distributed in IR 2 . The nodes are assumed to be slot-synchronized. In a randomly chosen slot, 
some nodes wish to transmit a packet. We assume that the locations x 1; x 2 , . . . of these transmitters 
follow an isotropic Poisson point process (PPP) $ := {x^}^ on R 2 with intensity^ X(x) being 
defined on IR 2 . Throughout this work we will refer by Xj to the random location of the z-th node 
as well as to the z-th node itself. 

Due to isotropy of $, X(x) is rotation-invariant [1_| and depends solely on the distance \x\ to 
the origin, i.e., X(x) = X(\x\e^ e ) for all 9 G [— 7r,7r)o For convenience, we define r := \x\. 

The following definition is a consequence of the fact that X(x) can be described as the resulting 
intensity after location-dependent thinning of a stationary PPP of some constant intensity lfT6l . 

Definition 1. The shape function F(r), reflecting the spatial shape of $, is defined by the relation 

\{x) = XF(\x\), (1) 

where F(r) > for all r > and max r {F(r)} = 1, and < A < oo is some intensity scaling 
constant. 

The restrictions are necessary to ensure that X(x) is non-negative and bounded by A every- 
where. The shape function F{r) is defined on IR> . 

We assume that each transmitter Xj has an intended receiver located at fixed distance d. In 
order to measure the (spatially-averaged) link performance in the network we define a reference 
link, cf. [4J; the reference link consists of a reference receiver placed at an arbitrary location 
yo G K 2 and of an associated reference transmitter placed at x , where x lies somewhere on 

3 Since the intensity of <J> directly relates to the average number of transmitters per unit area (at least throughout this work), 
we will use the terms "intensity" and "density" interchangeably without risk of confusion. 

4 In a very few situations, we will switch between Cartesian coordinates (x) and the corresponding polar coordinates (|x|e jfl ) 
when appropriate. We do not expect any confusion thereof. 
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Fig. 1. System model: The reference receiver is located at yo with distance \yo\ to the origin. The associated reference 
transmitter is located d units away from yo at location xq. Black dots represent interferers with random locations xi,Xa,. ... 
Shape function F(r) characterizes density of interferer set $ over distance r to the origin. 

the circle db(yo, d). Note that neither the receiver y nor the transmitter x are part of the point 
process $. By the Slivnyak-Mecke Theorem [1], the statistics of $ are not affected by the 
addition of the reference link. Note that the reference link x — > y is typical in the sense that 
it reflects the typical link performance for transmissions over distance d with receivers located 
at distance \y \ to the origin. 

We consider a path-loss plus block-fading channel with independent and identically distributed 
(i.i.d.) fading coefficients. The power path loss between two locations x,y G M 2 is given by the 
path-loss function £(\x — y\) := (c+ \x — y] )' 1 with path loss exponent a > 2. The parameter 
c > ensures boundedness of £ (modified path-loss model)@ The power fading coefficient 
between a transmitter at x and the reference receiver y is given by g x , where E [gj = 1 for all 
x. When appropriate we will drop the index x in g x . 

We assume that all nodes transmit with the same fixed transmit power and at a common 

5 Setting c = 0, we obtain the singular path-loss model, which is also often used in the literature. 
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information rate. The sum interference power at the reference receiver y is then given by 

l(2/o):=]Tgx^|x-?/o|). (2) 

It is further assumed that the sum interference signal is conditionally additive white Gaussian 
noise (AWGN), i.e., it follows a zero-mean Gaussian distribution with random variance I(t/ )- I n 
other words, conditioning on a spatial configuration of interferers Xi,x 2 , ■ ■ ■ as well as on each 
transmitter's channel gain to the reference receiver g XlV0 , g X2 y i • • •> me sum interference signal 
is zero-mean Gaussian distributed with (deterministic) variance I(yo). 

The instantaneous signal-to-interference-plus-noise ratio SINR at the reference receiver y is 
given by 

where rj is the average signal-to-noise ratio (SNR) resulting from white thermal noise in the 
receiver hardware. We assume that the nodes employ strong (capacity-achieving) channel coding, 
e.g., concatenated Turbo codes, such that the outage event is a steep function of the SINR. We 
focus on the case where practical constraints preclude explicit transmitter coordination as well as 
feedback of channel state information (CSI) to the transmitter. In such a scenario the transmitter 
is blind to the instantaneous realization of both the interference power \(y ) and the channel 
fading g XQ which may eventually lead to link outage. 

Definition 2. The outage probability of the reference link x — > y is given by 

q(y ) :=P(SINR(j/ ) </?), (4) 
where j3 is a modulation and coding specific threshold. 

Remark 1. In a more general scenario (whenever the Poisson assumption on $ is dropped), the 
Palm probability F lxo (■) should be used in dU) instead of P(-), since although the reference 
transmitter x does not contribute to the interference, it may influence the statistics of $. The 
probability measure P !x ° (•) characterizes the statistical law of the point process $U {x }, while 
removing the point x from the event ■ [24J. It can be seen as a conditional probability given a 
point x which, however, is excluded from the event in question. 
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A. Comments regarding some assumptions 

As for every model, its applicability has to be discussed and the underlying assumptions have 
to be justified. In this section we briefly comment on the assumptions made in this work. 

Despite its analytic tractability, the PPP assumption limits the degree of applicability of the 
model to realistic scenarios. This is because the Poisson property allows nodes to be arbitrarily 
close to each other with non-zero probability and furthermore precludes any explicit medium 
access coordination between the nodes (only Aloha-type medium access). Other models, such as 
the Matern hard-core point process, can in general better capture certain (physical) inhibition or 
coordination aspects in wireless networks, but at the cost of rapidly loosing analytic tractability. 
However, we note that there has been recent advances in the modeling of more sophisticated MAC 
strategies, such as carrier sense multiple access (CSMA) lfT6l . [|22|. Il23l and local frequency 
division multiple access (FDMA) ll26l . that use Poisson approximations. These methods were 
shown to work well in the high-reliability regime (at small target outage probability), which 
represents the more relevant regime from a practical viewpoint. Moreover, [21] recently showed 
that even for an arbitrary regular configuration of interfering nodes, e.g., hexagonal cell structure, 
the resulting interference field converges to a Poisson interference field as the variance of the 
shadow fading increases. For these reasons, we argue that the Poisson assumption still remains 
worthwhile. Besides, it offers the great possibility to reveal the basic effects of non-homogeneous 
spatial distributions and boundary effects in an analytic way. 

The fixed distance assumption, which can be seen as a target transmission distance that is 
dictated by the network protocol, is commonly accepted, see for example [3J. Furthermore, it has 
no effect on the outage probability analysis since the particular locations of the other receivers 
are irrelevant to the interference experienced by the reference receiver. Nonetheless, the fixed 
distance assumption will be dropped in Section IV-BI to account also for connectivity issues. 

It is well known that interference is appropriately modeled as AWGN when direct-sequence 
code division multiple access (DS-CDMA) or orthogonal frequency division multiplexing (OFDM) 
techniques are employed at the physical layer G71 . Whenever the AWGN assumption may not 
hold, our results can be seen as a lower performance bound, since one could generally do better 
by exploiting the interference signal at the decoding stage, i.e., by using interference-cancellation 
techniques J6]|. 
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B. Spatial shapes chosen for illustrations 

At this point we would like to introduce the four exemplary spatial shapes used for illustrations 
and numerical evaluations throughout this work. They are appropriately chosen in order to reflect 
common scenarios in wireless networks and to help increasing the reader's intuition about the 
result. They are depicted in Fig. [2] 

• Scenario a), finite network: This scenario reflects the basic property all practical networks 
share: the network area is finite, or equivalently, the node density tends to zero for sufficiently 
large distances to the origin. In order to allow comparability with prior works, it is assumed 
that the density first remains constant over a wide range as in the stationary case. At the 
network boundary the density then starts to decay rapidly until it becomes zero. 

• Scenario b), urban with hotspot: In urban scenarios there sometimes may exist small areas 
with very high traffic, i.e., communications hotspots. They are typically found in commercial 
areas or other public places, compromising many densely -and sometimes dynamically- 
deployed wireless architectures^ serving a large number of possibly mobile users. We model 
such a hotspot scenario by "adding a plateau of density" around the origin on top of an 
already existing level of density corresponding to the urban deployment. This level then 
decays to zero with increasing distance to the origin to reflect finiteness of the network. 

• Scenario c), scattered decentralized network: There are certain types of networking appli- 
cations that preclude a detailed network layout for the reason of hostile environments or 
limited geographic access. For instance, large sensor networks are sometimes created by 
airdrop which results in a highly scattered spatial distribution of devices. We model such a 
behavior by an exponentially decreasing density around the origin. 

• Scenario d), carrier sensing in decentralized networks: This scenario is found in decentral- 
ized networks with transmitters employing carrier sensing to avoid excessive interference; 
suppose a transmitter located in the origin o is granted access to the medium. Consequently, 
other potential transmitters directly surrounding this transmitter defer their transmission as 
they sense the medium as busy, while other potential transmitters farther away sense the 
medium as free and therefore start to transmit, cf. Section IV-Cl for more details. As a result, 
the density of active transmitters around the considered transmitter at o behaves effectively 

6 E.g., pico and femto cells in 3G/4G multi-tier networks, or Wi-Fi access points. 
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Fig. 2. The four exemplary spatial shapes considered in this work. Scenario a) (finite network), Scenario b) (urban with 
hotspot), Scenario c) (scattered decentralized network) and Scenario d) (carrier sensing in decentralized networks). 

as indicated by Scenario d) in Fig. [2l 

III. Interference Analysis 

We now study the interference statistics at the reference receiver y . The analysis will first focus 
on the case of general fading, i.e., with an arbitrary distribution. For the purpose of measuring 
the local throughput, the most valuable quantity of interest would be the distribution of SI NR. 
However, closed-form expressions exist only for a very limited number of channel models, see 
e-g-> 0- To overcome this problem, statistical techniques that bound the interference distribution 
have been proposed recently [3J. These techniques rely on the first moment of the interference 
as well as on the dominant interferer phenomenon, and are exploited in Section IIII-AI Here, we 
specifically make use of some characteristics of F(r). In Section lTlI-Bl we then concentrate on the 
case of exponentially distributed power gains (Rayleigh fading) and derive the Laplace transform 
of the interference distribution. The Laplace transform uniquely determines the distribution itself 
and is furthermore very useful for computing outage probabilities as we will see later on. 
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A. Arbitrary fading 

1) First moment of the interference: The first moment of the interference E [l(yo)] measured 
at the reference receiver y can in general be obtained by 

F(\(y )>z)dz. (5) 

Obviously, one would have to know the distribution of \(yo), which unfortunately is known in 
closed-form only for a few cases of stationary point processes [7]. A remedy to this problem is 
given by the Campbell Theorem flTJ , which allows us to derive the first result: 

Theorem 1. Let /(r) := dF(r)/dr[] c > 0, a = 2 and \y \ > 0. If lim F(r)r u < oo for some 

r— »oo 

v > 0, then 

E[\(y )] = \A 2 (y ,c)<oo, (6) 

where A 2 (y , c) is given by 

(I 1 2 POO 2 i I 1 2 \ 

F(0) asinh^f^ + / /(r) asinh 7 " ^"f ' dr . (7) 

A proof is given in Appendix El The condition lim F(r)r u < oo for some v > is necessary 

r— >oo 

for E[l(y )] to exist. The function A 2 (yo,c) in © has an interesting interpretation: it can be 
seen as the interference-driving function as it determines the interference up to a scaling factor. 
Besides, the first term in © can be interpreted as the interference field associated with the origin 
o, while the second term accounts for the particular spatial shape by effectively "adding up the 
interference in a differential way" according to f(r). 

Remark 2. If the reference receiver is located in the origin (\y \ = 0), the asinh-term in A 2 (yo, c) 
has to be replaced by log(r 2 + \y \ 2 + c), cf. Identity 2 in Appendix lAl 

Glancing at the second term in ©, we note: 

Corollary 1. When F(0) = 1 and /(r) < for all r e R + , F(r) can be interpreted as a 
complementary cumulative distribution function (CDF) with respect to a random distance r to 
the origin, yielding 

C -^° |2 ' - r asinh r2 + c 7 l5ol 2 1. (8) 



A 2 (yo, c) = -7T asinh— — + 7rE 
2 1 2/o I Vc 



2 1 2/o I \/c 



7 In cases where F(r) is not differentiable, e.g., F(r) = 1 — H(r) where H(r) is the Heaviside function |28 |, the differentiation 
is understood within the context of distribution theory. 
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Corollary Q] states that the integral in © can be seen as an averaging of the interference 
with respect to a random distance r. Such a representation may be appropriate when analyzing 
networks with a priori unknown or fast- varying spatial configurations, for which a CDF is then 
used to model their spatial shape. 

Corollary 2. Let < lim,,.^ F(r)r u < oo, for some v G (0, 2]. Then, the expected number of 
interferers E[|$|] = 2ir\ ' rF(r) dr = oo but the expected interference E [l(yo)] < °°- 

The intuition behind Corollary [2] is that, although the expected number of nodes contributing 
to the interference is unbounded, the network remains sufficiently sparse such that the first 
moment of the interference remains bounded. Note that for a PPP, if the expected number of 
interferers is unbounded, this implies that the number of interferers is a.s. infinite which can be 
verified by studying the Laplace transform of the PPP HI, ll29ll . Applying the Markov Inequality^ 
P( z > z) < ^E [z] [30J, this in turn means that the number of nodes is a.s. infinite while the 
interference remains a.s. finite. This particular finding is somewhat remarkable since it rearranges 
the commonly-accepted perception, stating that whenever a = 2 and the number of interferers 
is a.s. infinite, the interference is a.s. infinite as well JH, 0. This perception indeed holds for 
stationary point processes but does not hold in general for non- stationary point processes, as 
demonstrated by Corollary [2] 

Theorem 2. Let lim^oo F{r)r u > 0, for every v — > 0. Then, \(y ) = oo a.s. 

A proof is given in Appendix [Bj Theorem [2] shows that whenever F(r) decays at most 
logarithmically, the interference is a.s. infinite. In particular, this includes the stationary case 
(since lim^oo F(r) > 0) which is consistent with the literature [7|. Combining Corollary [2] and 
Theorem [2l we observe that for asymptotically decaying F(r) there exists a transition between 
sparse and dense networks. This transition determines whether or not the interference is a.s. 
finite in a non- stationary Poisson network with a.s. infinite number of interferers. 

Remark 3. By setting g = 1, the pure path loss case is also covered by the above results. 

We now treat the case a = 4. 

8 which shows that E[l(yo)] < oo implies that a.s. l(yo) < oo. 
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Theorem 3. Let /(r) := dF(r)/dr, c > and a = 4. Then, 

E[l(y )] = AA 4 (y ,c)<oo, (9) 

where Al(2/o, c) is given by 



a i \ 71 ( T?f \ f 2<H{«(r,c,2/o)} 
Ai(y , c) = — = F(r) arctan- 



r=0 



29 c t{/<r,c,2/o)} , . , m , 
/(r) arctan- — — dr (10) 

1 - \K{r, c,y )\ l 



%Vc \ 1 - \K{r,c,y )\ 2 

and K(r, c, y ) is given by (l34l) in Appendix |Al 

A proof is given in Appendix E The term At(yo, c) in (TTOl) can be again interpreted as the 
interference-driving function. 

Corollary 3. Let F(r) = 1. Then, /(r) = (stationary PPP). Then, by carefully taking the limits 

2^{ K (r,c )2 /o)} J" " 

lim arctan 



1 - \K(r,c,y )\ 2 



2, a = oo, 



we recover the well-known result E [\(yo)] = for the stationary PPP, cf. JT). 

2J Upper and lower bounds on the interference distribution: As mentioned in the beginning 
of Section Unl there already exist some techniques for bounding the tail of P (1(2/0) > 2). A 
simple upper bound can be constructed by applying the Markov Inequality and using the results 
from Theorem \T\ and Theorem [3l Since this procedure appears to be straightforward, we will 
first focus on the (more interesting) construction of a lower bound. First, we recall the definition 
of subharmonic functions. 

Definition 3. (Subharmonic functions OTl Ch. 2]): Let G C i" be an open set and let h(x) be 
a function twice continuously differentiable on G. If Ylk=i > on G, then h(x) is called 

k 

subharmonic on G. 

If F(r) is convex in a certain (one-dimensional) region, then the intensity \(x) is subharmonic 
on the corresponding (two-dimensional) region. Such a behavior may be often found at the 
network boundary, e.g., Scenario a) and b), or when the shape function exhibits a tail, e.g., 
Scenario c) and d). In this case a lower bound on the tail probability P (1(2/0) > z) can be 
derived: 
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Theorem 4. Let \(x) be subharmonic on G C IR 2 and let y G G. Denote by f(x) the maximum 
radius for which the closed ball b(x, t) is contained in G, i.e., f(x) = max{t G M + : 6(a;, t) C G}. 
Then, 

(/•^(s/o) \ 
-27i\F{\y \) j rF(g> z{c + r a )) drj . (11) 

A proof is given in Appendix El Note that subharmonicity includes the case of harmononicity. 
The construction of the lower bound in Theorem |4] basically builds on the so-called "dominant 
interferer" phenomenon introduced in flU, where it was also reported that the resulting bound is 
fairly tight. However, in our case the tightness of CCD) strongly depends on the second derivative 
of F(r) and may not be guaranteed. 

Corollary 4. Let z < -. For the pure path loss model (g = 1), CCD) reduces to 

P (Kvo) >z)>l-exp (-7rAF(|y |) min |f% ), g - c)S}) . (12) 

The restriction z < - is necessary to allow the closest interferer to be dominant, otherwise 
this bounding technique would yield the trivial lower bound F(\(y ) > z) > 0. 
Using the Markov Inequality 11301 , we obtain the simple upper bound 

P(lfeo) >z) < -A a (y ,c), (13) 

z 

for the cases a = 2 and a = 4, where A 2 (yo,c) and A 4 (?/o,c) are given by © and (flOl) . 
respectively. 

B. Rayleigh fading 

The case of Rayleigh fading channels is of special interest. It has been reported in several works 
01, CD that closed-form results for outage probability can be obtained via the Laplace transform 
of the interference distribution when the channel gains g follow an exponential distribution^ 
For this reason it is important to know the Laplace transform of the distribution of l(j/o), i.e., 
>C|( yo )(s) = E [e~ sl( - ? ' - ) ] . Due to the fact that the spatial distribution of nodes is non-stationary, it 

'Theoretically, it suffices that only the reference transmission xo —5- yo is subject to Rayleigh fading in order to exploit the 
Laplace transform of I (yo). It is however more appropriate to assume the same fading statistics for all channels due to symmetry 
reasons. 
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is very useful to derive C\( yo )(s) for general y and F(r) in order to obtain information about 
the interference also in the spatial dimension. 

Similarly to Section IIII-AI we focus on the two cases a = 2 and a = 4. While the latter 
scenario is often found in practice, the former one deserves a reasonable motivation: although 
being analytically tractable, the practical relevance of such a scenario (Rayleigh fading and 
a = 2) might seem questionable first; a = 2 is usually observed in free-space wave propagation 
environments while the Rayleigh fading model is appropriate for describing non-line-of-sight 
(NLOS) wave propagation with many reflected multipaths impinging at the receiver. It turns 
out that there truly exists urban NLOS scenarios with considerably small path loss exponents 
(a ~ 2.6) as demonstrated in 11321 . Depending on the geometry of objects in the proximity of 
the receiver, the received signal may therefore still undergo severe small-scale fading. Hence, 
the results for a = 2 may serve as a "theoretical limit" of what can be expected roughly in 
Rayleigh fading environments with small path loss exponents. 

Theorem 5. For exponentially distributed fading gains g (Rayleigh fading), the Laplace transform 
of \(y ) at y is given by 

£|(s/o)( s ) = exp(-\s A a (y , s + c)) , (14) 

for the cases a = 2 and a = 4, where v4 2 (?/o,c) and Al(?/o,c) are given by © and (flOl) . 
respectively. 

A proof is given in Appendix|Bl Note that for a = 2 and < lim^oo F{r)r u for every v — > 0, 
we have that C\( yo )(s) = for all s. This in turn implies \(y ) = oo a.s. which is consistent with 
Theorem [2l 

Remark 4. Setting F(r) — 1 for all r 6 M + and c = 0, we recover the well-known result for the 
homogeneous case with a = 4: C\( yo )(s) = exp(— A^v^)' c f- 0- 

IV. Outage Probability Characterization 

Outage probability as defined in © is an important physical layer metric for characterizing 
the performance of the network; clearly, if the outage probability is undesirably high the whole 
communication process will be degraded irrespective of the amount of engineering effort on the 
upper layers. For this reason outage probability is often considered as a bottleneck metric. It is 
thus worth to study the outage probability in detail. 
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Fig. 3. Outage probability 5(3/0) vs - \yo\- Parameters were d = 10, /3 = 0.5, c = 1, A = 0.001. The spatial shape function 
F(r) was chosen according to scenario a). Marks represent simulation results. 



A. Outage probability 

For arbitrary fading channels the bounds on the tail probability of \(y ) may be used to 
bound the outage probability, see for example [4]. In the sequel we will focus on the Rayleigh 
fading case and discuss the impact of the spatial shape on the resulting performance. As already 
mentioned, the Laplace transform of the interference power can be used to compute the outage 
probability. We therefore note: 

Corollary 5. The outage probability q(yo) at y for exponentially distributed fading gains g 
(Rayleigh fading)is given by 

q(y ) = l-C Kyo) (/3(c + d a ))e-n, (15) 

for the cases a = 2 and a = 4, where C\( yo ) is given by (fl4l) . 

The technique for obtaining Corollary \5\ is well-known in the literature, cf. 01, 0. It builds on 
the Laplace transform £,\r yo \(s) and on exploiting the fact that the g are exponentially distributed. 
For Rayleigh fading channels, we are now able by means of (fl"5l) to measure the exact outage 
probability at an arbitrary location y and for an arbitrary spatial shape function F(r) satisfying 
the given restrictions. 
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Fig. [3] shows q(yo) vs. \y \ for a = 2 and a = 4. The spatial shape function F(r) was 
chosen according to scenario a). In the noise-free case (r) = oo) the outage probability decreases 
monotonically with increasing distance to the origin. Furthermore, we observe that for a = 2 the 
outage probability q(yo) is higher and its slope is less steep than it is for a = 4. This is because 
for a = 2 the individual interference contributions decay more slowly over distance than they do 
for a = 4. As a result, the interference is no longer dominated by only a few nearby interferers 
but it is determined by the many nodes, including those relatively far away from y . 

When receiver noise is considered (r) = 10 dB) the behavior of q(yo) changes considerably: 
while q(yo) is on the same order as in the noise-free case around the center of the network, 
both curves converge to a constant outage probability level as the reference receiver entirely 
leaves the center of the network. In fact, in this boundary region outage is primarily due to 
bad fades as the density of interferers becomes negligible small, thus rendering the performance 
noise-limited rather than interference-limited. This transition -from the interference-limited to 
the noise-limited regime- can be precisely tracked owing to the model at hand; for example, 
Fig. [2] suggests that the noise-limited regime commences somewhere around \y \ m 500, while 
Fig. [3] reveals that this is not the case at least for a = 2 (\y \ > 800) for the reason explained 
above. 

B. Exact vs. approximate model 

Up to this point, it is not yet very clear how much can be gained by the interference model 
derived in this work. Clearly, the model at hand is precise and accurate since it properly captures 
boundary effects and other non-homogeneities, though at the cost of computing an additional 
integral due to the arbitrary shape function F(r). In order to quantitatively measure its benefits, 
we therefore compare our model to a much simpler one that uses a somehow natural way to 
approximate the non- stationary -property of the interference locally around the reference receiver 
y . In this simpler model, it is basically assumed that the interference field at y originates from a 
stationary PPP having constant intensity XF(\y \). This essentially means that the network- wide 
spatial distribution of interferers is approximated locally by the density at location y . 

For the stationary case the outage probability for Rayleigh fading is well-known [7j. With the 
above approximation the approximate outage probability q(yo) at y is then given by 

q(y Q ) = 1 - exp (-AF(|y |)d 2 /3^ csc ^ . (16) 
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The intuition behind (fl6l) is that the exact outage probability is approximated by taking the 
outage probability expression corresponding to the stationary case and modulating the intensity 
A by the spatial shape function F(r). In order to measure the difference between q(yo) and q(yo) 
we define the following metric. 

Definition 4 (Log-divergence). The log-divergence (LD) is defined as 



7(y„):=A-Mogl^4M. (17) 

1 - 9(2/0) 

The LD 7(1/0) quantifies the ratio of the exact and approximate success probability^ on a 
logarithmic scale for arbitrary y . The normalization by A is necessary to remove the dependency 
on A so to measure the divergence resulting from the spatial shape only. Since the actual 
divergence technically scales with A, 7(1/0) should be interpreted as the logarithmic success 
probability divergence per density A. For large positive 7(2/0) me approximation overestimates 
the true outage probability, while for large negative 7(2/0) outage probability is underestimated. 
The approximation is accurate whenever 1 7(2/0) I is small. Using the result from Corollary |5] we 
are now able to evaluate the degree of accuracy of our model for a specific case. 

l0 The success probability is defined as 1 — q(yo)- 
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Corollary 6. Let the g be exponentially distributed (Rayleigh fading) and let c = 0. The LD for 
the case a = 4 is then given by 



Fig. |4a] shows the LD 7(2/0) vs. \yo\ for the shape functions F(r) introduced in Section 
III-Bl It can be seen that the LD exhibits an oscillatory behavior that depends on the degree 
of variability of F(r). For example, F{r) in scenario a) and c) changes comparably slow, and 
thus the corresponding LD shows only weak oscillations. The LD for scenario a) remains low 
(7(2/0) ~ 0.80) over a wide range which suggests that the local approximation works well in this 
case. As for scenario c), however, the LD is large around the origin which is due to the fact that 
the outage probability is highly overestimated as the (exponential) decay of F{r) to the right- 
hand side is neglected by the approximation. A similar effect can be observed for scenario d): 
the outage probability is highly underestimated around the origin because the increasing density 
to the right-hand side is neglected. The LD for scenario b) exhibits rich oscillations due to a 
stronger varying shape function. These oscillations are strong particularly in the transition region 
(10 1 < r < 10 3 ) as illustrated by the inner graph of Fig. Hal The latter also holds for scenario 
a) since the LD increases by approximately 3.8 times around the transition point r = 10 2 in this 
case. 

In order to understand how outage probability is actually affected by this divergence the 
relative approximation error 



is illustrated in Fig. [4b] for the example A = 1CT 3 . The relative approximation error is an important 
measure for quantifying the outage probability deviation, and thus for characterizing the expected 
QoS variability. Fig. [4b] underlines the observations made in Fig. [4a] the approximation works 
well for scenario a) while for the other three scenarios the approximation is relatively loose. 
Especially for the two scenarios b) and d) the relative approximation error S(y ) can be very 
high over a wide range: for scenario c) the relative error is between 1% and 10% around the 
hotspot. For scenario b) the relative error is approximately 10% around the density mid-level. In 
case of scenario d) the approximation completely fails around the origin. Interestingly, for the 
first three scenarios the relative approximation error starts to increase rapidly when passing the 




(18) 



<%o) 



Q(yo) - q(vo)\ 
q(yo) 



(19) 



November 21, 2012 



DRAFT 



22 



transition point r = 10 2 . In contrast, the relative approximation error decays for scenario d) as 
\y \ becomes large which is a result of F(r) being asymptotically constant, cf. Fig. [2} Note that 
the zeros in Fig. [4b] correspond to the zeros in Fig. |4al 

V. Applications of the Model 

We now discuss how our model can be applied to problems in wireless networks with a non- 
stationary spatial node distribution. The question of how to measure local throughput in such 
non- stationary networks is addressed first. 

A. Local transmission capacity 

As discussed in Section II- Al the original notion of TC can no longer be applied in the case 
of non-stationary networks. It is however possible to study the throughput of such networks in 
a way similar to using the TC metric. This in fact requires a refinement of the TC definition 
which is given in the following: 

Definition 5 (Local transmission capacity). The local TC is defined as the maximal density 
of concurrent transmissions in an infinitesimal region around the point x subject to an outage 
probability constraint e, i.e., 

c(x,e) := A(x,e)(l-e). (20) 

The TC and its "local" counterpart have similar meaning, except that the latter is location- 
dependent: For a given location-dependent density X(x) and target outage probability e, c(x, e) 
yields the TC in an infinitesimal region around the point x. In this way the local TC implicitly 
takes into account the spatial shape of the node distribution. For isotropic node distributions the 
local TC only depends on \x\, and hence it gives the TC in an infinitesimal region around all 
points that have distance \x\ to the origin. For Rayleigh fading, the local TC c(x, e) is obtained 
by solving (fl"5l) for A. Otherwise, the upper and lower bounds on the interference distribution 
can be used to obtain bounds on the local TC. Note that the perhaps most desirable feature of 
the TC, which is the ability to benchmark different transmission protocols, is preserved by the 
local TC. 

In order to demonstrate the use of the local TC we revisit the famous question whether to 
employ narrowband or wideband signaling in interference-limited networks, cf. [|5]|. 
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Fig. 5. Local transmission capacity ratio C DS , Q ' e : vs. lyol and spreading gain M. Parameters are: /3 = 1, a = 2, d = 10, 
c = 0, Rayleigh fading model. The scenario b) was chosen. 



Example 1 (Frequency hopping (FH)-CDMA vs. direct sequence (DS)-CDMA in decentralized 
Poisson networks). In H it has been reported that for the stationary PPP model, the TC gain 
of FH-CDMA compared to DS-CDMA scales as M 1_ q, where M is the processing gain. The 
intuition behind this result is that, as the path loss exponent increases, it gets more and more 
preferable to avoid interference than to average it out. However, for very small path loss exponents 
(a — > 2) this result suggests that the benefit of FH-CDMA vanishes irrespective of the processing 
gain M, i.e., the gain is M 1_1 = 1. In fact, it turns out that the correct scaling behavior for 
a = 2 cannot be properly captured by the analysis in 01 due to the underlying stationarity 
assumption. Using the model presented here, we are able to correct the above result for a = 2 
using the local TC metric. This is demonstrated in Corollary [7J for the case of Rayleigh fading 
and y = o, i.e., with the reference receiver located in the origin. 

Corollary 7. Let a = 2, c = and r] = oo\^j For exponentially distributed fading gains g 
(Rayleigh fading), the local TC gain of FH-CDMA over DS-CDMA at o is given by 

c Db (o,e) A 2 {o,(3d 2 ) 

"The requirements c = and r] — oo are not imperative for the derivation of this Corollary. However, they were chosen for 
the reason of comparison with |4l . 
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A proof is given in Appendix [Bj Surprisingly, this result shows that FH-CDMA is still better 
than DS-CDMA for a = 2 as opposed to the result in [4J; the local throughput gain is always 
greater one and scales with logM. Although Corollary [7] assumes the reference receiver to be 
located in the origin, we conjecture that this scaling behavior is not changed when the reference 
receiver is located at an arbitrary location. Fig. [5] supports our conjecture. 

The plot also demonstrates a fundamental property of the local TC: the performances of 
different schemes can be analyzed and compared over the network area. It provides additional 
information about the health of the network in space domain; this potentially may help net- 
work operators and designers in identifying performance bottlenecks and to allow for proper 
engineering, e.g., to balance QoS among all nodes irrespective of their location in the network. 

B. Spatial throughput 

Besides the TC, spatial throughput is another metric that is widely used in the field of 
network analysis, e.g., [[16]. In the stationary case, the spatial throughput is usually given as 
the average number of successful transmissions per unit area. Here, the normalization to unit 
area is appropriate since -by the stationarity property- the local performance is statistically 
the same everywhere in the network. This normalization, which essentially gives the spatial 
throughput the character of a density, is furthermore necessary to ensure that this throughput 
metric is bounded; clearly, the absolute throughput is a.s. infinite since infinitely many nodes 
contribute to it. A key advantage of this metric is that it is analytically tractable. It can be 
written as the density of transmissions times the (Palm) probability of success associated with 
this density. 

As discussed earlier, real networks are usually finite, and hence non- stationary which results in 
non-homogeneous statistics and boundary effects. These effects unfortunately cannot be captured 
by the spatial throughput metric since it requires stationarity. A more suitable metric is proposed 
in the following: 

Definition 6 (Sum spatial throughput (SST)). The SST is defined as the expected number of 
successful simultaneous transmissions, i.e., 




(22) 
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Note that a sufficient condition for the SST to be bounded is that the (expected) number of 
transmitters is a.s. finite. However, it may be theoretically possible that infinitely many nodes 
transmit simultaneously but only a finite number of these transmissions is successful. 

Similar to the spatial throughput metric, the expectation operator inside the SST is computed 
using the Campbell Theorem JT), lfl6ll . l|29l . In contrast to the local TC, which provides spatial 
information about the local throughput, the SST computes the sum local throughput taking into 
account the spatial non-homogeneity. In particular, boundary effects are properly captured as 
opposed to the traditional spatial throughput metric for the stationary case. 

The success function 1{ X successful}, indicating whether transmitter x is successful, can be chosen 
arbitrarily to include additional outage-inducing effects, e.g., energy-limitations, dis -connectivity 
or secrecy outage. The definition of the SST given by ((22)) is rather general. In the following, we 
show how the isotropy property can be exploited to compute the SST for a scenario of special 
interest. 

Theorem 6. With the underlying network model and success function l{siNR(y)>/3}> the SST p 
becomes 



A proof is given in Appendix El 

In the following, we will give an optimization example in order to demonstrate how the SST 
can be used to solve problems in network design. The chosen example does by no means cover 
all aspects and impairments of a practical network and is meant to serve as an illustration. 
However, the fixed distance assumption initially made in Section |n] will be dropped to account 
also for connectivity issues. The considered optimization example can be easily extended to 
more sophisticated problems. 

Example 2 (Optimization of SINR target (3). We consider a wireless network with scattered nodes, 
e.g., a sensor network created by airdrop. In one arbitrary snapshot, the set of all nodes can be 
divided into the set of transmitters and the set of potential receivers. Without any modification, we 
assume that the transmitter locations are characterized by the isotropic PPP $ with density X(x) = 
XF(\x\). Define by {yi}^ the set of potential receivers and let the yi,y2, ... be independently 
distributed on R 2 with average density A r (r) = A r F(r), where A r > 0. We assume slotted 
Aloha as the MAC, rendering $ and {yi}^x independent (bipolar model), cf. ifToll . In such a 




(23) 
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scenario connectivity at distance d is no longer guaranteed for every node. We further assume 
that the routing protocol employs a nearest neighbor strategy, i.e., among all potential receivers 
every transmitter chooses its receiver such that d is minimized. Since the distance between every 
transmitter-receiver pair is now random, we write d. For independently distributed points, the 
distribution of d between a given point and its nearest neighbor is well-known, cf. ifToll . 

As for the optimization, we would like to know the optimal modulation and coding parameters 
such that the Shannon throughput is maximized under the SST metric. This corresponds to 
maximizing the product plog 2 (l + (3) over the target SINR 0. Thereby, a transmission from a 
transmitter x to its nearest receiver y over a distance d = |x — y| is considered successful if 
the receiver y is not in outage, i.e., l{siNR(x, y )>/3}LJ We have altered the notation SINR(y) — > 
SINR(x, y) to indicate the random distance d between transmitters and receivers. The optimal 0* 
is then obtained by solving 



0* = argmax <^ log 2 (l + /3)E l{ S iNR(x, y )>/3} 
Applying the Campbell Theorem 01, lfl6l . E9l . we rewrite the expectation in (l24"l) as 



P 



Li 2 



E [l{SINR(x,y)>/3}] A(x) dx. 



(24) 



(25) 



Next, we condition the nearest receiver of transmitter x to be located at y = y which is equivalent 
to conditioning on d = d and = 9, where 6 is the random angle enclosing transmitter x and 
its receiver y. This is easily seen by noting that y = x + de^ e in polar coordinates. Hence, 



P 



P(SINR(x,y) > 0) 



1-9(2/) 



X(x) dx. 



(26) 



Conditioned on the the distance d, the term P (SINR(x, y) > 0) now depends only on the 
interference \{y) as well as on the channel between transmitter x and its nearest receiver 
y, since the particular location x no longer matters. It can be computed using the methods 
presented in Section |IV] to obtain the (conditional) outage probability 1 — q(y). To evaluate 
the inner expectation we would have to know the distribution of 1 — q(y) with respect to 



It may happen that two transmitters select the same node as their nearest receiver. In such a case, we assume that the 
receiver tries to decode both signals independently of each other, i.e., by treating the other signal as pure noise even though it 
may have been decoded already. One the one hand, this assumption is necessary to ensure analytic tractability, on the other it 
is tolerable since we are not investigating achievable rates. 
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TABLE I 

The three cases used for the optimization example. The parameter 7710 denotes the average SNR at 

DISTANCE d = 10. 





A 


Ar 


rjio 


Q 


Case I 


10~ 3 


10 -3 


12 dB 


2 


Case II 


10 -3 


10~ 2 


12 dB 


2 


Case III 


10~ 2 


10 -2 


12 dB 


2 



random angle 0, which has a non-uniform distribution in the non-stationary case. To get rid 
of this analytic intractability, we approximate the outage probability at y averaged over by 
the outage probability experienced at x. The intuition behind this approximation is that on 0- 
average, the distance of receiver y to the origin is in the order of the distance of transmitter x 
to the origin. Since the outage probability depends on the distance of the receiver to the origin, 
this approximation will hence reflect the 0-averaged outage probability at y up to acceptable 
accuracy. Finally, we obtain the following simplified expression 

/*oo /»oo 

p~2n\ rF(r) (1 - q(r, d)) dP(d = d) dr, (27) 
Jo Jo 

where we altered the notation q(r) — > q(r, d) to highlight the dependency of the outage probability 
on d. 

Fig. [6] shows the plog 2 (l + 0) vs. (3 for different variations of the average SNR, density of A 
and A r as shown in Table H First of all, we observe a tight matching of theoretical and simulation 
results, indicating that the approximation made above is indeed justified. The plot furthermore 
reveals some interesting implications concerning the right choice of 0: for example, the f3* 
maximizing plog 2 (l + 0) is more sensitive to the ratio j- and less sensitive to the average SNR, 
even when taking into the account the noise-limitedness at the network boundary. Furthermore, 
the effect of a limited SNR is insignificant when the density of transmitters is high which is 
consistent with our intuition, though this evidence is now quantitatively confirmed owing to the 
model at hand. 

Alternatively, the optimization problem can be reformulated as a spatial-aware optimization 
problem, i.e., after node deployment by optimization over all functions /3(r). 
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/3 dB 

Fig. 6. Maximization of p\og 2 (l + 0) for the Rayleigh fading case and scenario c), e.g., decentralized sensor network. The 
system parameters are shown in Table U Solid lines represent the no-noise case (7710 = 00), while dashed lines represent the 
case with average SNR 7710 = 12 dB at distance d = 10. Marks represent simulation results. 



C. Interference in CSMA networks 

Besides slotted Aloha, other MAC protocols, such as CSMA or local FDMA, represent effec- 
tive techniques for reducing excessive interference generated by nearby nodes. These techniques 
coordinate transmissions locally to inhibit nodes from accessing resources that are already in 
use. In order to study decentralized networks with such inhibition mechanisms while ensuring 
analytic tractability, powerful methods based on non-homogeneous Poisson approximation have 
been proposed recently ifToll . [|22|. Il26l . When such protocols are transmitter-initiated, e.g., 
transmitter sensing for CSMA or transmitter orthogonalization for FDMA, the resulting spatial 
distribution of interferers becomes non-homogeneous but remains isotropic around the inhibiting 
transmitter. This is because potential transmitters around the inhibiting transmitter are kept silent 
while others located farther away are likely to transmit. In contrast, the interference field at the 
associated receiver is not isotropic at all in general, as illustrated in Fig. [7J In order to maintain 
analytic tractability the interference field at the receiver is however often assumed isotropic as 
if the receiver were co-located with the inhibiting transmitter ll22l . Il26l . It turns out that the 
accuracy loss due to this simplification strongly depends on the distance between the inhibiting 
transmitter and the associated receiver; clearly, as the receiver gets closer to the inhibiting 
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Fig. 7. Snapshot: resulting interferer set after CSMA mechanism with transmitter-sensing. Grey dots represent inhibited 
interferers and black dots represent active interferers. An interferer is inhibited whenever it senses an ongoing transmission, here 
transmission of xq. Interference field around receiver yo is not isotropic. 



transmitter this simplification becomes more and more justified. However, we are unable to 
quantify this accuracy loss as a function of distance unless time-consuming experiments and/or 
simulations are conducted. 

Fortunately, our model is capable of evaluating this loss since it perfectly fits into such an 
isotropic and non-homogeneous situation. Before demonstrating that, we briefly recall the basic 
idea of this particular modeling technique for CSMA networks for ease of comprehension. The 
reader is referred to jl22j. 11261 for further details. 

Example 3. Let the set of potential interferers be initially distributed according to a stationary 
PPP of density A. Assume that the reference transmitter x is located in the origin and the 
associated reference receiver y is located d units away from xq. The inhibition mechanism can 
be appropriately embedded into the model using a two-step approach [|22|: First, the large-scale 
density of active interferers is derived using a Matern-type 2 point process model [1], thereby 
capturing the inhibition effect on a "macroscopic" level. The large-scale density \ e is then given 
by 

1 — e ^ a' 

Xe = — ^ 2T^^' (28) 

7rr(l + J)A~- 

where A > is the CSMA power sensing threshold. As an intermediate step, we condition on 



November 21, 2012 



DRAFT 



30 



the fact that x is granted access to the channel and we are now interested in the statistical 
properties of the interferers surrounding x after this conditioning. Clearly, the conditioning 
affects the activity status of the potential interferers as the transmission of x may be sensed by 
the potential interferers. To overcome this analytic intractability the set of potential interferers is 
approximated by an non-homogeneous PPP in the second step. Here, the problem of correctly 
capturing the interactions between the nodes surrounding xq and xq itself is virtually transformed 
into a location-dependent thinning of the PPP. The resulting small-scale density A s around the 
reference transmitter x is then given by 

A,(r) = A,(l-exp(-Ar Q )), (29) 

where the term 1— exp (— Ar a ) can be seen as the probability that an interferer at distance r to xq 
does not sense the ongoing transmission of Xq. The density A s (r) is monotonically increasing in 
r since it is more likely that far away potential interferers will be unable to sense the transmission 
of xq. The behavior of A s (r) can be appropriately described by the spatial shape function of 
scenario d) with density scaling factor \i, cf. Fig. [2l The sensing threshold effectively controls 
the trade-off between interference reduction on the one hand and spatial reuse on the other. 

Finally, the density A s (r) in (1291) is used to describe the interference around the reference 
receiver y , although it reflects the interference at the reference transmitter xq. This simplification 
procedure increases analytic tractability at the cost of losing accuracy in later computations. The 
amount of this accuracy loss is demonstrated in the following for the case of Rayleigh fading 
channels, a = 4 and c = 0. 

We use the same notion in (fl~9~l) to measure the relative accuracy loso as a function of the 
distance d, i.e., 

5(d) = , \y \ = d, (30) 

QW, d) 

where the second argument in q now highlights the dependence of the outage probability on 
d. Fig. [8a] shows the loss in accuracy vs. d for different 0. As for the inhibition threshold, a 



13 We use the term accuracy loss rather than approximation error in order to avoid misunderstandings: the two-step approach 
explained above for embedding the CSMA mechanism into the model already inherits an error due to the Poisson approximation. 
However, we are interested in the additional loss in accuracy that is induced by the simplification procedure described above. 
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(a) Relative accuracy loss (b) Absolute accuracy loss 

Fig. 8. (a) Relative accuracy loss 8(d) vs. transmission distance d for different /9. Parameters are a = 4, c = 0, A = 10~ 3 
and A = —50 dB. (b) Absolute accuracy loss \q(o) — q(yo)\ vs. d for different /3 and A. Parameters are a = 4, c = 0. A was 
chosen such that the spatial throughput is maximized under the maximum outage probability constraint e = 0.25. 



moderate value of A = —50 dB was choseno The density of potential transmitters is A = 1CT 3 . 
It can be seen that depending on the distance d, large relative errors occur. Unfortunately, they 
are spread over a wide range of d with their highest value concentrated around d = 10; a 
distance that is often used in the modeling of decentralized networks 0), H|. Interestingly, the 
relative accuracy loss decreases as the target SINR (3 increases. This can be explained by the 
dominant interferer effect: the region of dominant interferers around both x an d yo increases 
with (3 which in turn increases their overlapping, resulting in a higher outage correlation. In other 
words, the sum interference created by nodes farther apart becomes more and more relevant as 
increases. In this case the interference field is no longer dominated by a few nearby nodes but it 
is dominated by many far away nodes that have a small relative difference between the path loss 
to Xq and to y . On the opposite, at low (3, e.g., when employing DS-CDMA, the interference 
is dominated by a few nearby interferers, thus rendering the interference experienced at xo and 
yo less correlated. 

14 A larger A would increase spatial reuse at the cost of increasing outage probability, while a smaller A would decrease 
outage probability at the cost of decreasing spatial reuse. 
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Similar to (1301) , now the absolute accuracy loss \q(o)—q(y )\, where \y \ = d, is shown in Fig. 
[8b1 vs. d. In order to ensure practical relevance as well as to avoid certain combinations of system 
parameters that would lead to pathological cases, e.g., very high outage probability, we chose the 
parameter A as follows: for each d choose the A that maximizes the spatial throughput subject 
to an outage constraint of e = 0.25. This simple optimization problem can be easily solved 
numerically ll33l . The absolute accuracy loss is of special interest since it is necessary to track 
the absolute deviation from targeted outage probability, i.e., the QoS requirement. The plot clearly 
demonstrates what was already discussed above: the accuracy loss increases with decreasing (3 
(error > 20% for (3 = 0.1). Lowering the density of potential interferers A effectively shifts high 
accuracy losses towards larger d. Note that assuming the reference receiver y to be located in 
the origin yields a too optimistic outage probability. In this example, this consequently means 
that in the dense regime, where q(yo) = e (horizontal lines), the actual outage probability is 
q(yo) — q(o) times higher than the targeted value of e = 0.25 (+21.6% for (3 = 0.1, +9.6% for 
13 = 1 and +1.68% for (3 = 10). 

In conclusion, we note that if the network is driven at high f3, i.e., spectral efficiency of 
transmission is high, the accuracy loss induced by assuming that the reference receiver sees 
exactly the same interference field as the reference transmitter may be negligible. Unless the 
transmission distance d is ultimately small, i.e., short-range transmission, the exact location of 
the reference receiver y relative to the reference transmitter x must be considered in all other 
cases in order to avoid high losses in accuracy. 

VI. Conclusion 

We extended prior work on the interference modeling and throughput analysis of wireless 
Poisson networks by assuming an isotropic but not necessarily stationary spatial distribution of 
nodes. In such non- stationary networks, which are usually found in practice, the performance 
is location-dependent and the spatial distribution of nodes may have an arbitrary shape. In this 
work, we analyzed the interference statistics and the outage probability as a function of (i) an 
arbitrary location inside the network and of (ii) an arbitrary shape of the node distribution. For 
the path loss exponents a = 2 and a = 4 we obtained closed-form expressions which led to 
some interesting observations. For example, for a = 2 it is possible that the interference is 
almost surely finite while the number of contributing interferers is almost surely infinite; thereby 
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demonstrating the existence of a sparse and a dense regime. This particular result sharpens prior 
statements on the interference characterization for the case a = 2. The interference analysis 
was conducted for arbitrary fading models (including the pure path loss case) as well as for the 
Rayleigh fading model. Moreover, we presented a lower and upper bound on the tail probability 
of the interference for arbitrary fading, where the lower bound is not limited to a = 2 and a = 4. 

Furthermore, we proposed two performance metrics, namely the local transmission capacity 
and the sum spatial throughput, which are suitable for measuring the local throughput in such 
non- stationary networks. The local transmission capacity is a refinement of the well-known 
transmission capacity, which fails to characterize throughput in networks with non- stationary 
spatial node distributions. This refinement allows for comparing and benchmarking different 
protocols over the network area, and provides spatial information about the local throughput. 
The sum spatial throughput counts the average number of successful transmissions in the network, 
thereby accounting for boundary effects and other non-homogeneities in the spatial distribution 
of nodes. The sum spatial throughput may be of special interest for optimizing finite networks, 
which due to their non- stationary nature, cannot be properly described using the traditional spatial 
throughput metric. 

We applied our interference model developed in this work to these two metrics and discussed 
their importance with respect to some problems in network design. For example, using the 
local transmission capacity we were able to show that for very small path loss exponents 
(a — > 2) the throughput gain of FH-CDMA over DS-CDMA scales logarithmically with the 
processing gain, and thus is non-vanishing as opposed to prior investigations. Our model recovers 
accuracy losses coming from local approximations that neglect boundary effects and/or non- 
homogeneous deployments in existing networks. Seeing the model as a design tool, it may be 
of significant importance for the correct dimensioning of system parameters in order to satisfy 
QoS requirements. We also showed that the developed model can also help to better describe 
networks that employ CSMA sensing at the transmitter. 
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Appendix A 
Integral Identities 



Identity 1. If a > \b\, a,b el, 



d0 



vrP„ 



Va 2 -b 2 



l (a + bcos^ (a?-b 2 )^ ' (31) 
where P n (x) is the n th -Legendre polynomial. We will be using the O th -Legendre polynomial 
given by Pq(x) = 1. 

Identity 2. Let ai,a 2 ,a 3 G R, where a 3 > 0. Define R :— ai + a 2 t 2 + a 3 t 4 , T = 4a!a 3 — a 2 ,. 
Using the substitution t — >■ t 2 , we have 

T > 



Vai + a 2 t 2 + a 3 t 4 



dt = < 



asinh 203^02 _ 



(32) 



log(a 3 t 2 + si 



2 /' 



T = 0. 



Identity 3. Let ai, a 2 G R, where ai > 0. Then, 

2t 

cti + if 2 + a| ~ 2ta2 cos 0) 2 
where 




7r 29Uk;(£, ai, (Z2)} 

atan 



2^01 1 — ai, 02)! 



2 • 



n(t, ai, a 2 ) := 



(33) 



(34) 



V / (v / ^ + j(i 2 + ^)) 2 + 4t 2 a 2 
Identity 3 deserves a proof: 

Proof: The basic idea is to decompose the integrand into partial fractions and to apply 
Identity [T] and |2l yielding (l33l) after some algebraic manipulations. Note that according to P4ll . 
(13TI) and (|32l) hold only for real- valued parameters. However, we verified that they also hold for 
complex- valued parameters. □ 



Appendix B 
Proofs 



A. Proof of Theorem [7] 

We want to compute the expectation 

E 



5>^(ix-2/o|) 



(35) 
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where the expectation is with respect to all interferer locations x G $ and all channel gains g x . 
Since the expectation operator linear, we can compute the expectation with respect to all g x first, 
i.e., E [g x ] = 1 Vx G $. Applying the Campbell Theorem, yields 

E[l(y )] = A / £(\x-y \)F(\x\)dx (36) 

Changing to polar coordinates and exploiting the isotropy property, we can rewrite (1361) as 

r°° r 2r F(r) 
E [!(!/„)] = A/ / , / o i i lo o i i ^ d <f> dr - 07) 



o jo 



c + (r 2 + 1 2/0 1 — 2 r \yo\ cos 

We then apply Identity Q] to the inner integral of (1371) to obtain 

r°° 2r F(r) 

mVo)] = ^ 7( =. 2 / 2 ; = 2 , dr. (38) 

Jo V( c + r + \Vo\) +4r 2 |y | 2 

Finally applying Identity [2] to d38l) and verifying the convergence of the upper limit using the 
constraint lim F{r)r u < oo, for some v > 0, yields the result. □ 



B. Proof of Theorem \2\ 

To prove that \(y ) is a.s. infinite, we analyze its Laplace transform E [exp(— s\(y ))} and 
check if it is zero for all s. In the PPP case, the Laplace transform of the interference field is 
given by 0, [[El 



E[exp(-sl(j/ ))] = exp -E 



f h _ e -«g<(l»-i»l)) x(x) dx 



(39) 



For ([39l to become zero, the integral must not converge. We write 

f (i _ e -*i\-vo\)\ X (x) dx > I *et(\* -Vo\) X{x)d ( 40) 

where (a) follows from the inequality < 1 — e~ z for z > — 1 [28 J. Inserting the path loss 
function £(\x — y \) = (c + \x — j/o| 2 ) -1 into the right-hand side of (|40l) yields 

At the upper limit of the inner integral the integrand behaves as F(r)/r. So, the condition 
< lim^oo F(r)r u < oo, where v — \ 0, is sufficient for the divergence of the integral. Because 
the Laplace transform of \(y ) becomes zero in this case, this concludes the proof. □ 
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Fig. 9. Schematic illustration of the proof of Theorem [4] Shaded area around yo represents the largest disc contained in the 
subharmonic region G. 



C. Proof of Theorem \3\ 

The proof is analogous to the proof of Theorem [T] except for the integration part. Here, the 
integral Identity [3] is used instead. We further exploit the fact that max r {F(r)} = 1 to ensure 
convergence of the integrals. □ 

D. Proof of Theorem |?] 

We follow the idea of dominant interferers which was introduced in [4] to bound the tail 
probability ¥(\(y Q ) > z): let A{z) := {(x,g x ) G $ x R + : gj(\x - y \) > z) denote the set of 
all interferers, where each one taken by itself already results in the event \(y ) > z. Clearly, 

lfeo)= gx^(|x-2/ |) + gx^(|x-yo|), (42) 

(x, gx )eA(z) (x, gx )e.4(z) 

where the set A is complimentary to A on the same domain, i.e., A(z) := {(x, g x ) G $ x M + : 
gx^(|x — I) < z }- Note that both sums in (|42|) are non-negative. Hence, we can now write 

> z) > p I Yl gx £ (i x - f°D ^ 2 

\(x,gx)e^l(z) 
= P(l^)l>0) 

( = ) l-exp(-A(^))), (43) 



where (a) follows from removing the "non-dominant" part in (1421) . (b) follows from the definition 
of A(z) and (c) is a consequence of |^4(-z)| being Poisson distributed with mean A(A(z)). Using 
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lfT6l Corollary 2.1.2], A.(A(z)) in (g3) can be computed as 

A(A(z))=[ F(g> z{c+\x -y \ a ))X(x)dx. (44) 

We now translate $ by the vector y to obtain a ?/o-centric coordinate system, cf. Fig. [9j After 
switching to polar coordinates we obtain 

A(A(z))= r¥ (g > z (c + r a )) X(y + re j<t> ) d0dr. (45) 
Jo Jo 

We now exploit the fact that X(x) is subharmonic around y in a region G. Let f(y ) denote the 
maximal radius for which b(y , f(y )) is contained in G. Then, (|43T) can be bounded as 

(a) MoT) /-27T 

A(.A(z)) > / rF(g> z(c + r a )) X(y + re 1 *) d<pdr 
Jo Jo 

(b) fr{Vo) 

> 27rA(y ) / rP(g > z(c + r a )) dr 



(c) r (yo) 

= 2nXF(\y \) r F (g > z (c + r a )) dr, (46) 



where (a) follows from limiting the upper integration limit to f(y ). Inequality (b) is a conse- 
quence of subharmonicity ll35l Ch.10], and (c) follows from the CQ) in Definition 1. Inserting 
(1461) in (O yields the result. □ 



E. Proof of Theorem \5\ 
We write 

exp(-sg x £(|x- y \)) 

= exp (- J ^ ^1 - £ g - 3/oD) J A(ar) dx 

where (a) follows from algebraic manipulations and the i.i.d. property of the g x . (b) follows from 
the probability generating functional and the Laplace functional of a PPP IfToll . After noting that 
C g (s) = (1 + s)^ 1 for exponentially distributed g, the integral is computed using Identity \T\ and 
Identity [2] for a = 2 and Identity [3] for a = 4. 

Note that (a) in the proof holds for general point processes and some approximation techniques 
for computing the right-hand side already exist |f36l| . The (b) part is for PPPs only. □ 



C 



l(2/o)l 



(a) 



E, 



ri E g 
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F. Proof of Corollary [7] 

Solving (TT3T ) for A and multiplying by 1 — q(y ) yields the local transmission capacity 

-lQg(l-6)(l-6) 

f3d 2 A 2 {y ,(3d I ) 

after substituting q(yo) — > e. We are interested in the case yo = o. Using Identity 2 (case T = 0) 
we find that 



POO 

A 2 (o, fid 2 ) = -F(0) log (3d 2 - / f(r) log(r 2 + /3d 2 ) dr. 

Jo 



(48) 



Assuming that nodes employ pseudo-noise sequences, FH effectively thins out the point process 
of interferers by M (interference avoidance, A/M), while DS reduces the interference by a factor 
of M (interference averaging, 0/M), cf. [4J. Hence, using (|47T) and (|48T) the ratio ^ DS |°'^ can be 
written as 

c FH (M = M ^ 2 , ±2) 

c DS (o,e) (3d 2 A 2 (o, (3d 2 ) M n ' M J 

_ ^2 (0, |^ 2 ) 
A 2 (0,/3d 2 ) 

( „) -F(0) log |^ 2 _ /^/(Qlog^ + y) dr 

A 2 (o,/3d 2 ) A 2 (o,/3ci 2 ) 
(m , F(0)logM 1 f°° r 2 + (3d 2 , 

A 2 (o,,5d 2 ) A 2 (o,/3d 2 )J r 2 + j^d 2 

where (a) follows from (fl"5l) and Remark [2] and (b) follows from algebraic manipulations. Now 

we show that the integral in (|49l) is 0(1). Assuming \ f(r)\ < oo (which is reasonable in practical 

networks), we note that the integrand has no singular values. For M = 1 the logarithm becomes 

zero and so does the integral. Since log r + ^ < log(l + ^f-) for all M > 0, we therefore 

analyze the convergence of the integral 

rOO rOO 

J /(r)log(l + 2£)dr <max{\f(r)\}J log(l + ^) dr < oo. (50) 
Hence, the integral in (|49l) is bounded by a constant. □ 
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G. Proof of Theorem [6| 
We write 



E 



^ ^ l{x successful} 
.xe<3> 



(a} 

(J 
(J 
(e) 



E successful}] A(x) dx 



e$ [1{sinr(v)>0>] dp (y = i/k) d ?/ -M^) dx 

E$ [l{siNR( y )>/3}] dP (y = y\x) \{x) dx^j dy 
P$ (SINR(y) > 13) dP (y = y|ar) X(x) dx) dy 



(y = y\x) X(x) dx dy 



1 - g(y)) 



;i - q(y))Hy) d v- 



(51) 



(a) is due to the Campbell Theorem, (b) is obtained by noting that a transmitter x is successful if 
the intended receiver at y is not in outage. From Section UH we know that y is placed by random 
translation of x according to some probability kernel dP (y = y\x). (c) follows from Tonelli's 
Theorem 03 and (d) follows from E [l{xeA}] = P (X G A), (e) follows from CES) and the fact 
that q(y) is independent of x. □ 
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